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The  spectra  of  Nd^"*"  and  Er^'*'  in  CaWOtT  are  reanalyzed  by  using 
free-ion  wave  functions  to  determine  new  phenomenological  crystal 
field  parameters,  Ratios  of  [bj^mKEr ) ,'bjj^  (Nd))  are  shown 

to  be  approximately  independent  of  k and  m. where 

^km  “ ^km  [ ^ ^ “ '^k  )]  This  result  reinforces  the  previous 

assumption  that  derived  may  be  factored  according  to 
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ion  independent.  Certain  assumptions  regarding  excited  radial 
wave  functions  in  addition  to  these  results  lead  to  theoretical 
T and  <r>Cr^>^j  n'^'^'n'  values  for  nf  = 5d,  5g  for  all  the  lanthan- 
ides. Together 'with  <r*^>HF  O]^  values  obtained  from  the  liter- 
ature, one  may  derive  ■3'^d  calculate  theoretical  electric 

dipole  oscillator  strengths  for  any  lanthanide-host  combination 
where  the  Aj^rn  “^re  determined  by  a lattice  sum  over  the  host’s  con- 
stituent ions.  Tables  required  for  performing  these  calculations 
are  given. 


A/cr 


1 . INTRODUCTION 


In  the  first  three  reports  of  this  series. 


method  was 


established  for  predicting  the  crystal  field  parameters,  Bj— , of  triply 


ionized  lanthanides  in  crystals.  These  Br 


for  even-k  define  the 


crystal  field  that  determines  the  observed  Stark  splittings."  a’ 
rationale  was  given  for  factoring  the  of  the  lanthanide  ions  in  a 
given  host  according  to 


^km  ^ ^ ^k)\m  ' 


where  the  are  coefficients  in  the  spherical  tensor  decomposition  of 
the  electrostatic  field  at  the  impurity  ion  site  and  are  specifically 
impurity  ion  independent.  Thus,  it  was  argued,  one  can  calculate  the 
A^  by  performing  a lattice  summation  over  the  constituent  ions  of  a 
given  host  material  and  determine  the  B,  for  any  lanthanide  according 
to  equation  (1) . If  variations  from  host  to  host  in  the  radial 
integrals  and  shielding  factors^  can  be  assumed  small,  then, 

presumably,  B,  can  be  derived  for  any  impurity  ion-host  system 

/k\ 

once  <^r  y and  are  known  for  the  lanthanides  in  some  convenient  host. 


Phenomenological  B 


for  even-k  determined  by  the  effective 


spin-orbit  Hamiltonian  (ESOH)  procedure^  for  six  lanthanides  in 


^N.  Karayianis  and  C.  A.  Morrison,  Rare  Earth  Ion-Host  Lattice 
Interactions,  1.  Point  Charge  Lattice  Sum  in  Scheelites , Harry  Diamond 
Laboratories  TR-1648  (October  1973) . 

2/V.  Karayianis  and  C.  A.  Morrison,  Rare  Earth  Ion-Host  Crystal 
Interactions , 2,  Local  Distortion  and  Other  Effects  in  Reconciling 

Lattice  Sums  and  Phenomenological  B , Harry  Diamond  Laboratories 
TR-1682  (January  1975)  . 

P.  Leavitt,  C.  A.  Morrison,  and  D,  E.  Wortman,  Rare  Earth 
Ion-Host  Crystal  Interactions , 3.  Three-Parameter  Theory  of  Crystal 

Fields,  Harry  Diamond  Laboratories  TR-1673  (June  1975)  . 

J.  Freeman  and  R.  E.  Watson,  Phys.  Rev.,  127  (15  September  1962), 

2058. 

^P.  Erdos  and  J.  H.  Kang,  Phys.  Rev.  B,  ^ (1  November  1972),  3393. 

^N.  Karayianis,  J.  Chem.  Phys.,  £3  (15  September  1970),  2460. 


CaWOi.  were  used^  to  test 
simplest,  nontrivial  forms 


and 


so  that  derived  (as  opposed  to  phenomenological)  B,  (t;qri)  fot  even-k 

are  functions  of  the  three  parameters,  i,  q,  and  n-  The  t parameter  in 

this  theory  corrects  for  inadequacies  in  Hartree-Fock  radial  wave 

functions,  and  the  values  of  t which  were  found  to  be  about  0.75  are 

consistent  with  the  t necessary  to  adjust  Pr^"*"  Hartree-Fock  values  for 

Slater  integrals  to  empirical  values.^  The  q parameter  is  the  effective 

charge  on  the  anion  in  the  covalent  radical  of  the  scheelite-type 

crystals  considered  (for  example,  the  effective  oxygen  charge  in  the 

wo 4 radical  in  CaW04) , and  n is  the  ratio  of  the  effective  distance  of 

that  charge  from  the  cation  to  the  anion-cation  nuclear  distance  as 

determined  by  x-ray  or  neutron  diffraction.  This  so-called 

three-parameter  theory  of  crystal  fields  was  implemented^  to  determine 

how  well  one  can  fit  derived  to  phenomenological  B,  for  the  six 

km 

lanthanides  Nd,  Tb,  Dy,  Ho,  Er,  and  Tm  in  CaW04 . It  is  important  to 

establish  how  well  the  electrostatic  field  of  the  host  material  can  be 

represented  by  the  A,  (qn) , because,  if  the  A,  for  even-k  are  adequate 

fcm  Kin 

to  predict  the  B,  , then  the  A,  for  odd-k — important  in  intensity 
km  km 

calculations — also  may  be  presumed  adequate.  It  is  presently 

impractical  to  obtain  the  latter  by  any  other  method. 


Karayianis  and  C.  A.  Morrison,  Rare  Earth  Ion-Host  Crystal 
Interactions,  2.  Local  Distortion  and  Other  Effects  in  Reconciling 
Lattice  Sums  and  Phenomenological  B,,  Harry  Diamond  Laboratories 
TR-1682  (January  1975) . 

^R.  P.  Leavitt,  C.  A.  Morrison,  and  D.  E.  Wortman,  Rare  Earth 
Ion-Host  Crystal  Interactions , 3.  Three-Parameter  Theory  of  Crystal 
Fields,  Harry  Diamond  Laboratories  TR-1673  (June  1975)  . 
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this  premise.  It  was  concluded  that  the 
for  the  and  A^^^  in  equation  (1)  are 


c\  -k  / k\ 

/ /hf 


(2) 


A = A,  (qn)  , 
)on  km  ^ 


(3) 


1 


. .t 


r 


Point  charge  lattice 


sums 


were  performed,  and  t , 


and  1 were  varied  to  fit  derived  to  phenomenological  B,  independently 

for  the  six  ions  in  CaW04.  It  was  then  shown,  by  fixing  the  effective 

oxygen  charge  to  q = -1.09  and  its  fractional  distance  from  tungsten 

to  h = 0.977,  that  a best  f for  each  of  the  six  lanthanides  could  be 

determined  that  gives  substantially  the  same  fit  to  phenomenological  B, 

km 

T,  q,  and  n were  varied  independently.  The 


as  previously  obtained  when 
fact  that 

(-1.09,  0.977)  represents  a qualified  success  for  the  theory,  but  the 


fits  to  the  B,  of  each  ion  could  be  made  with  the  same  A, 

km  km 


rms  deviations  remaining  between  derived  and  phenomenological  B,  are 

km 

still  too  large,  ranging  from  13  (for  Er)  to  135  cm“'  (for  Tb) , to  give 
one  confidence  in  the  predicting  capabilities  of  the  theory. 

p.  ^ 

In  this  report,  the  latest  values  for  quantities  necessary  to  derive 

the 

discrepancies . 


Bj^  are  used  in  an  attempt  to  isolate  the  source  of  these  large 


new  analysis  using  free-ion  wave  functions  is 
performed  on  the  spectra  of  Nd  and  Er  in  CaWO(7  to  determine  more 

iai< ) 

accurate  (©^factors  for  these  ions.  |>.A  linear  interpolation  to  all  the 
lanthanides  of  these  t,  along  with  smoothed  values  for  Freeman  and 
• and  linearly  interpolated  values  for  the  a of  Erdos 


Watson ' s 


HF 


and  Kang,^  are  used  to  determine  where 


P,^(T) 


= - %)  ' 


(4) 


so  that 


(5) 


P.  Leavitt,  C.  A.  Morrison,  and  D.  E.  Wortman,  Rare  Earth 
Ion-Host  Crystal  Interactions,  3.  Three-Parameter  Theory  of  Crystal 
Fields,  Harry  Diamond  Laboratories  TR-1673  (June  1975)  . 

‘*A.  J.  Freeman  and  R.  E.  Watson,  Phys.  Rev.,  127  (15  September  1962), 
2058. 

^P.  Erdos  and  J.  H.  Kang,  Phys.  Rev.  B,  £ (1  November  1972),  3393. 


Values  are  derived  also  for  quantities  necessary  for  calculating 

transition  probabilities  such  as  ',r^^  for  p = 1,  3,  5,  and  7 

and  A , = E „ - E,„  for  nil  = 5d  and  5g. 
ni  nH  4f  ’ 

The  results  given  in  this  report  set  the  stage  for  calculating 
spectra  and  electric  dipole  transition  probabilities  effectively  from 
first  principles  for  all  the  lanthanides.  These  calculations  will  be 
performed  for  many  of  the  current  experimental  host  materials  to 
establish  an  initial  competency  in  predicting  potential  laser  properties 
of  new  materials.  Results  of  this  project  will  be  the  topic  of 
subsequent  reports  in  this  series.  While  it  may  be  overly  optimistic  to 
expect  that  one  can  accurately  predict  spectra  once  the  crystal 
coordinates  of  the  host  material  are  known,  it  is  hoped  that  some 
guidance  can  be  given  as  to  which  materials  show  greatest  promise  as 
laser  materials. 


2.  DERIVING  CRYSTAL  FIELD  PARAMETERS 


The  first  step  in  deriving  the  p for  the  lanthanides  according  to 

^ 

equation  (4)  requires  a determination  of  \r  / and  a . We  have 

Hr  K 

obtained  fifth-degree  polynomial  expressions  in  N that  fit  the  values 

of  \^r  / calculated  by  Freeman  and  Watson^  for  Ce,  Pr,  Nd,  Sm,  Dy,  Er, 
Hr 

and  Yb  to  better  than  0.24  percent  accuracy.  Letting 


/' 

\ 


k 


r 


-L 

mss  r 


km 


(6) 


where  N is  the 
lanthanide,  one 


number  of  f electrons  for  a given  triply  ionized 
gets  values  for  the  r^^^  given  in  table  I.  These 


“♦a.  J.  Freeman  and  R.  E.  Watson,  Phys,  Rev.,  127  (15  September  1962), 
2058. 
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TABLE  I . 


VALUES  FOR  r,  IN  UNITS  OF  A 
km 


WHERE 


HF  ■ ^m'’km'^ 


^ 2m 


4m 


6m 


0 3.79955(-l) 

1 -5.I  395K-2) 

2 8.3209K-3) 

3 -8.92996(-k) 

L 5.023021-5) 

5 -1.12879(-C) 


3.k5k3l (-1) 
-9.0M55(-2) 
t .78303(-2) 
-2.05603(-3) 
I .2058L(-M 
-2.78753(-6) 


6.56k96(-l) 
-2.38090(-1) 
5. 129k5(-2) 

-6.15291 (-3) 
3.69091 (-9) 

-8.6L990(-6) 


.VOTE;  Sumbers  in  parentheses  are  powers  of  10. 


formulas  give  values  for  Gd (N  = 7) 

of  (r*^)  = 0.223,  0.124,  and 

- k 

0.148  in  units  of  A for  k = 2,  4, 
and  6,  respectively,  whereas 
Freeman  and  Watson  obtain  0.220, 
0.119,  and  0.138  by  their  Hartree- 
Fock  method.*  The  differences  be- 
tween corresponding  values  indi- 
cate the  degree  of  accuracy  that 
one  may  expect  in  extending  the 
polynomial  fits  to  the  lanthanides 


for  which  radial  integrals  were  not  obtained  by  Freeman  and  Watson.  For 

/ k' 

completeness,  table  II  gives  the  \r  , for  all  the  triply  ionized 

nr 

lanthanides  calculated  from  equation  (6)  by  using  the  rj^^  values  in 
table  I. 

TABLE  II.  VALUES  FOR  t,  (^'^)  AND  ‘^'^HF  ' ‘’k^ 

DERIVED  CRYSTAL  FIELD  PARAMETERS  ARE  GIVEN  BY  B,  = P,  A,  FOR  GROUND 

km  k km 

CONFIGURATIONS  OF  TRIPLY  IONIZED  LANTHANIDES 


fon 

S 1 

1 ,,  ^ ^ ' 

Ce 

1 

0.7693 

0.3360 

0.2709 

0.4659 

0.6757 

0.0254 

-0.0422 

0.1841 

0.7536 

2.3417 

Pr 

2 

0.7597 

0.301*1 

0.2213 

0.3459 

0.6667 

0.0272 

-0.0421 

0.1756 

0.646** 

1 .875** 

Nd 

3 

0.7500 

0.2803 

0.  1882 

0.2715 

0.657S 

0,0290 

-0.0420 

0.1706 

0.5776 

1-5897 

k 

0.71*0} 

0.2621 

0. 1655 

0.2247 

0.6488 

0.0308 

-0.0418 

0. 1679 

D.5339 

1 .4213 

Sn 

5 

0.7306 

0.2472 

0. 1488 

0.1929 

0-6398 

0.0327 

-0.0417 

0.1668 

0.504? 

1 ’710 

Eu 

6 

0.7210 

0.231*7 

0.1353 

0.1686 

0.6309 

0.0345 

-0.0415 

0. 1666 

0..836 

1-25C3 

Gd 

7 

0.7113 

0.2232 

6.1237 

0.1477 

0,6220 

0.0363 

-0.0414 

0.  1668 

0.4656 

1.1873 

Tb 

e 

0.7016 

0.2129 

0.1131 

0. 1287 

O.6I3O 

0.0381 

-0.0413 

0.1673 

0.4??0 

1.1232 

Oy 

9 

0.69'9 

0.2033 

0,1037 

0.’11« 

0.6041 

0.0399 

-0.041  1 

0. 1681 

0.M34I 

1 .06 '4 

HO 

to 

0.6S23 

0.191*5 

0.0951* 

0.0981 

0.5951 

0.0418 

-0.0410 

0.1692 

0.4217 

1 .0119 

Er 

1 1 

0.6726 

0.1865 

0.0883 

0.087'. 

0.5861 

0.0436 

-o.otoe 

0.1706 

0.4126 

0.9826 

Tn 

12 

0.6629 

0.1790 

0.0820 

0.0787 

0.5772 

0.0454 

-0.0407 

0.1722 

0.4053 

0-964? 

Yb 

13 

0.6532 

: 0.1717 

0.0753 

0.0681 

0.5683 

0.0472 

-0.0406 

0.1737 

0.3938 

0.9120 

Values  for  the  were  obtained  by  linear  interpolation  by  using  the 

Erdos  and  Kang^  values  for  Pr  and  Tm.  One  gets 

*Multiply  the  values  of  Freeman  and  Watson,  which  are  in  atomic 
units,  by  (0 .529172)^  to  obtain  units,  where  the  angstrom  unit 
A = 10"®  cm. 

^P.  Erdos  and  J.  H.  Kang,  Phys.  Rev.  B,  £ (1  November  1972),  3393. 
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2^ 


. .1. 


w 


02 


0.6846  - 0.00895N 


(7a) 


04  = 0.02356  + 0.00182N 


(7b) 


06  = -0.04238  + 0.00014N  , 


(7c) 


which  generate  the  values  given  in  table  II. 

Values  for  the  t given  in  table  II  are  calculated  by  the  formula 

T = 0.75(1.0387  - 0.0129N)  , (8) 

which  was  obtained  by  a linear  interpolation  of  values  obtained  for 
Nd^'*'  and  Er^'*'  in  CaW04  as  follows:  Free-ion  basis  wave  functions'  for 
the  lowest  14  and  10  multiplets,  respectively,  for  Nd  and  Er  were 
obtained  by  using  the  "free-ion"  parameters  of  Carnall,  Fields,  and 


Rajnak . 


Crystal  field  subspaces  for  S4  symmetry  were  then  set  up  and 


diagonalized  for  the  crystal  field  Hamiltonian 


H = I,  B”  E.C,  (r.) 
X km  Ian  1 km  1 


where  the  C are  unnormalized  spherical  harmonics.  By  varying  the  B, 
km  km 

and  the  multiplet  centroids,  fits  of  4.56  and  3.80  cm“^  rms  were 

obtained,®  respectively,  to  51  (out  of  64)  and  38  (out  of  58) 

experimentally  identified  Stark  levels  of  Nd  and  Er.  These  parameters 

are  given  in  table  III,  columns  2 and  4.  If  the  phenomenological  B, 

km 

are  factorable  according  to  equations  (4)  and  (5)  where  the  A,  are 

km 

impurity  ion  independent,  then  b,  , defined  by 

km 


b,  = B, 
km  km 


(9) 


Chew.  Phys.,  49 


T.  Carnall,  P,  R.  Fields,  and  K.  Rajnak,  , 

(1968),  4424. 

E.  Wortman,  C.  A.  Morrison  and  N.  Karayianis , Rare  Earth  Ion-Host 
Lattice  Interactions , 5.  Lanthanides  in  CaWO^^,  Harry  Diamond 

Laboratories  TR-1794  (June  1977). 
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by  using  the 
different  ions 


two 


r and  values  from  table  II,  and  evaluated  for 

, will  be  related  according  to 

fb,  (Er)/b^  (Nd)V^’'  = T_/T^  , 

L km  km  J Nd  Er 

independently  of  k and  m.  The  b^^  and  values  given  in 

table  III  demonstrate  the  approximate  constancy  in  the  latter  values 
with  the  possible  exception  of  the  value  calculated  from  b20-  Giving 
double  weight  to  the  4,4  and  6,4  components,  we  obtain  approximately 

Er  Nd 

from  the  average  of  the  values  in  table  III,  column  6. 


TABLE  III. 


PHENOMENOLOGICAL.  B,  FOR  Nd3+  AND  Er^"^  IN  CaWO^  AND  CALCULATED 
. km 


Nd 

~ tr  Nd  'd<-f  vi*dl 

I 

6 

>kfn 

B t> 

Kp  kfn 

Nd^’E 

' B, ,, 0.252, 

-1.15.  0.9621 

■B 

20 

50J 

5244 

133  5610 

1 .0543 

393 

1 1 0 

40 

-878 

-4806 

-655  -7756 

1 .1271 

- 1062 

I8>< 

44 

1042 

5705 

809  9580 

1.1385 

1119 

77 

60 

-5.84 

-13  57 

-2.56  -28.14 

1-1293 

5 58 

q 

d64 

905 

3t99 

546  60J2 

I . M06 

734 

171 

i6m 

238 

841 

156  1715 

1 ,1261 

41 

197 

' 115 

wt  ighted 

1 4 1 f '*'4 

VJTL': 

Appr^pr  :j('  '.’til 

jt's  jn-  jsf^ 

!ron  tjbii-  tl.  pfrivf 

i flir-  ■ h 

fj‘i  jf  JiJtr  Ji ; f «.r*'nv'<*s  f.  twf.-n  drr: 

V''d  jnd 

ph»Tiom<‘nc..it)gic*i/  vaiues  for  .vJ  ifi  ijst  tor)  •.'Jiinns. 


To  obtain  values  for  each  t , we  performed  a lattice  sum  over  the 
CaW04  ion  setting,  in  units  of  proton  charge,  = +2  and  = -2  - 4q 
and  varied  the  oxygen  charge,  q,  its  effective  fractional  distance  from 


11 


tungsten,  g,  and  to  obtain  a best  fit  to  the  phenomenological  of 
Nd.  The  resulting  values  for  t,  q,  and  n,  which  are 

T^d=0-75,  (12a) 

q = -1.150  , (12b) 

n = 0.962  , (12c) 

gave  a fit  of  141  cm“^  rms  between  derived  (table  III,  column  7)  and 

0 + 

phenomenological  (table  III,  column  2)  for  Nd  in  CaWOu.  This 

"fit"  is  worse  than  that  obtained  to  the  phenomenological  B,  resulting 

Kin 

from  the  inherently  less-accurate  ESOH  method.^  Assuming  a linear 
variation  for  t over  the  lanthanides,  equations  (12a)  and  (11)  give 
equation  (8)  and  the  calculated  values  for  t in  table  II,  column  3. 
According  to  equation  (4) , one  may  then  calculate  the  values  for 
given  in  the  last  three  columns  of  table  II. 

The  main  ingredients  for  deriving  B,  for  any  lanthanide-host 

km 

combination  have  thus  been  obtained,  provided  that  one  assumes  that 

the  T,  approximately  host  independent.  The  final 

task  in  deriving  Bj^  in  a given  host,  under  these  assumptions,  is  to 

develop  an  adequate  model  of  the  particular  host  material  for 

calculating  the  A,  — the  q,h  model  assumed  in  the  three-parameter  theory 
km 

is  but  an  initial  effort  in  this  direction.  The  large  discrepancy  of 

141  cm”^  rms  between  derived  and  phenomenological  B,  for  Nd^"*"  in 

km 

CaW04  suggests  that  further  theoretical  work  is  required  to  determine 
the  ultimate  validity  of  this  theoretical  approach. 


3f?.  P.  Leavitt,  C.  A.  Morrison,  and  D.  E.  Wortman,  Rare  Earth 
Ion-Host  Crystal  Interactions , 3.  Three-Parameter  Theory  of  Crystal 


Fields,  Harry  Diamond  Laboratories  TR-1673  (June  1975) 
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3.  CALCULATING  ELECTRIC  DIPOLE  INTENSITIES 


Q 

To  calculate  electric  dipole  intensities,  one  must  know  A,  values 
, Km 

V 

for  odd-k,  r expectation  values  between  the  ground  4f  and  excited 

nS,  configurations,  where  I + 3 = odd,  and  energy  differences  between  the 

configurations  = E^j^  - Specifically,  the  matrix  element  M^(ij) 

that  determines  the  electric  dipole  intensity  of  polarization  u between 

N 

the  states  i and  j in  the  ground  f configuration  is  given  by 

= E V,<^>4f,n/^^)4f,ni^n^;M''m^ih^ 

Jkmnio 

X j [3][aJ  [j]  1 (3(0)1(0)U(0))  <Jl(0)k(0)  |3(0))  (13) 

{ W(3k31;iJ)<k(m)l(y) |j(m+u))  , 


where 

is  a unit  spherical  censor  for  N equivalent  electrons, 

Qc]  E 2x  + 1, 

(mj)  j2  (m2)  I j (m)^  is  a Clebsch-Gordan  coefficient,^^ 
and  W(abcd;ef)  is  a Racah  coefficient.^^ 

The  oscillator  strength  of  the  transition  is  given  by^'^^ 

P = X [4Trm  v,./ftlM^(ij)  , (14) 

y ei^Jy 


R.  Judd,  Phys.  Rev.,  127  (1  August  1962),  750.  Note  that  our  A'j^ 
are  Judd’s 

W.  Nielson  and  G.  F.  Foster,  Spectroscopic  Coefficients  for  the 
p^ , d^  and  f^  Configurations,  The  MIT  Press,  Cambridge,  MA  (1963). 

^^M,  Rotenberg , R.  Bivins,  N.  Metropolis,  and  J.  K.  Wooten,  The  3-j 
and  6-j  Symbols,  The  Technology  Press , Cambridge , MA  (1959) . See  their 
p.  1,  eg  (1.1)  to  relate  Clebsch-Gordan  coefficients  to  3-j  symlools  and 
p.  13,  eg.  (2.1)  to  relate  Racah  coefficients  to  6-j  symbols. 

^^M.  J.  Weber,  T.  E.  Varitimos  and  B.  H.  Matsinger,  Phys.  Rev.  B,  8_  (1 
July  1973)  , 47. 


2 2 

where,  in  terms  of  the  refractive  index  n,  X = n(n  + 2)  /9,  m is  the 

e 

electron  mass,  is  the  frequency  of  the  radiation,  and  -ft  is  Planck's 

constant  divided  by  2n . 

As  is  customary,  we  assume  that  mixing  between  the  4f  and  the 
excited  5d  and  5g  configurations  is  primarily  responsible  for  electric 
dipole  transitions.  To  obtain  approximate  values  for  the  radial 

integral  between  these  configurations,  we  assume  orthogonal i zed 
hydrogen-like  radial  wave  functions.  The  wave  functions  are  defined  by 
comparing  expectation  values  with  experimental  and  theoretical  data  from 
the  literature  wherever  possible. 

We  choose  4f,  5d,  and  5g  normalized  wave  functions  as  follows: 

= W^^r‘*exp  (-yr)  , (15a) 

= W^,(r^  + Ar**  + Br^)exp(-xr)  , (15b) 

bd  5d 

= Wj.  (r^  + Cr*^ ) exp (-xr)  , (15c) 

5g  5g 


where 


Jr  dr  ill  .4j  . - 6 . . 

n 1 3 ^3 


From  equation  (16) , the  normalizations  and  the  coefficients  A,  B,  and  C 
are  defined  in  terms  of  y and  x (see  app  A) . 

To  determine  y as  a function  of  N,  the  number  of  electrons  in  the 
N 

ground  f configuration  of  a triply  ionized  lanthanide,  ratios 

of  (r^).,  = 39/y^  were  adjusted  to  fit  similar  ratios 

calculated  by  Freeman  and  Watson  using  their  Hartree-Fock  wave 
functions.*^  The  quadratic  relationship 

y(N)  = 5.53358  + 0.228834N  - 0.00482594n2  a” ^ (17) 

J.  Freeman  and  R.  E.  Watson,  Phys . Rev.,  127  (15  September  1962)  , 

2058. 


gives  fits  of  1.03  percent  or  better  to  the  eight  calculated  ratios  of 
Freeman  and  Watson,  as  shown  in  table  IV. 


TABLE  IV.  RATIO  OF  FREEMAN  AND  WATSON  HARTREE-FOCK  EXPECTATION 
VALUES  FOR  SEVERAL  TRIPLY  IONIZED  LANTHANIDE  GROUND 
fN  CONFIGURATIONS 


N 

''hf 

y(NI 

39[v  N ] - 

(A-) 

(A‘M 

(A-  - ) 

lA-  ) 

1 

1 .201 36 

5.69850 

5.75759 

1 . 1 7698 

2 

1 .08628 

5.99263 

5.97199 

1 .09359 

3 

o 

o 

O 

6.22999 

6. 1 7665 

1 .02225 

5 

0.903F6 

6.57009 

6.55710 

0.90707 

7 

0.8I2S6 

6.92778 

6.89895 

0.81991 

9 

0.76992 

7.19052 

7.202IQ 

0.75186 

1 1 

0.70666 

7.92873 

7 .96662 

0.69951 

13 

0.65160 

7.73695 

7.69289 

0.65901 

\'OTE: 

Kt  t fCtivv  hydrotft'nic  u-jv*.* 

t unction 

iKiiii-nt  in 

— i-x;, 

lud-iFu  tic 

■ bt.'St  t' 1 1 Q(^)  = ‘ 

O.UHS  M.V 

- 0.  • 

, un  i >rti- 

s;r)n<iin>j 

rjtios  assLimin  ; ^ fx/)/' 

1 . 

To  determine  x,  calculated  ratios  of  ' r^) ^ ^ ^ were  fit  to 

4f , 5d  ■ 4f  ,5d 

similar  ratios  calculated  from  Grossgut's  wave  functions  for  Ce  , 

“Mi 

Pr^  , and  Nd  . A linear  variation  of  x according  to 

x(N)  = 2.48605  + 0.044175N  A" ^ (18) 

gives  a fit  to  Grossgut's  ratios  of  2.00  percent  or  better  as  shown  in 
table  V. 


The  complete  set  of  calculated  y,  x,  and  expectation  values  for  all 
the  triply  ionized  lanthanides  is  given  in  table  VI  where  columns  b,  c, 
d,  and  e are  normalized  to  column  a for  each  N. 

Grossgut,  Analytical  Wave  Functions  and  Oscillator  Strengths  for 
Nd^'^  and  Cu^'*' , ph.D  Dissertation,  Texas  Christian  University 
(1971)  (University  Microfilms,  Ann  Arbor  MI,  No.  72-7621) . 
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TABLE  V.  RATIO  OF  CROSSCUT'S  HARTREE-FOCK  EXPECTATION  VALUES  FOR 
Ce,  Pr,  AND  Nd  CROUND  CONFI CURATI ONS 


N 

(A2) 

'hf 

«(N) 

(r>) 

) 

1 

2.20785 

2.55597 

2.53023 

2.23898 

2 

Z.igiAB 

2.52291 

2.57ALO 

2.13591 

3 

2.0IA22 

2.6AA32 

2.61858 

2.03955 

NOT£: 

"iffective"  hydroqenic 

wave-function 

exponent  in  ‘ 

exp  <-XfifS) 

jssuminv  - exp  (•’y(N)r)  with  y(N)  from  table  IV,  linear  best  fit  x(S)  » 

2.'t8605  * 0,044175  N,  and  ratio  of  expectation  values  calculated  with  x(N) 
for  comparison  with  Grossgut' s values. 


TABLE  VI  . 

BEST-FIT 

EXPONENTIAL 

FAdORS 

IN  ~ exp 

[-y  (N)  r] 

AND 

pa 

V ’N‘ 

^'■Ni 

.1 

b 

C 

;A*  ' 

u--) 

(A-) 

(m  ) 

-I 

^ i 

1 

b.lblbS 

2.53023 

0.1)491 

2.23898 

n. 202621 

U. 766971 

3.2828’ 

2 

5.3715^ 

2.57^*^ 

9.60786E-2 

2.I35M 

0.18852? 

0.67IC99 

2.70:3 

} 

6.17665 

2.61858 

8.03255E-7 

2.0^955 

0.1 77361 

''.5961?.’ 

2.26565 

(« 

6. 371  7 

2.662  75 

6.9390IE-2 

1 .9512 

0. 168i.58 

0.536615 

1 .93739 

6.5571 

2 . 70693 

6.O636U-2 

1 .36996 

0.  I6I347 

O.'.oSSoS 

1.67-45 

h 

6.73285 

2.75H 

5.35826E-2 

1 .7953'* 

0.155688 

0.450039 

1 .47102 

7 

6-89895 

2.79528 

'..73655E-2 

1 .72685 

0.151232 

0.1.13399 

1 .30892 

7.05539 

2.S39'.5 

i..32078e-2 

1 .66396 

0.147794 

0.392475 

1.17853 

5 

7.20218 

2.80363 

3.93978E-2 

1 .60621 

0, 145236 

0.371216 

1 .07288 

10 

7.33933 

2.9278 

3.6273IE-2 

1.55314 

0.143456 

a. 3538s 

0.936835 

1 1 

7. '4668? 

2.97197 

3.37077E-2 

1 .50434 

0.I'.2373 

0. 339716 

0.916565 

12 

7.58.65 

3.01615 

3.)6041£-2 

1 .H^SkU 

0. 1419^47 

0.328^15 

0.359'93 

7.6928'. 

3.O6O32 

0.029886 

1 .‘.1808 

0.142127 

0.319578 

0.91255- 

- . .|f;j  , -X  , 

'.ir:  .*  r.jpi  t.ibl(ys  IV  arui  .in./  . 

V.llcuJ-ltl  '/  < 

..xpCK-rnr  i on  v'.i .'  i'’.*? 

j.'ifl  r.if.o.*?  ‘'roTr?  c.5t*.so  w.iv< 

f it.ct  !on.< . 

• r r-  ••>)!  ; 

V2  ..  17, 

" t"  •<!  "1 

, ■ <■><•■> 

<,> 

Vi  in-  . 1 

■■  >,. 

,•/ 

1-  . ;/)  ■- >I 

<'  >., 

h CO. 

The  final  factors  required  for  electric  dipole  intensity 

calculations  are  the  energy  differences  A , = E „ - E,,  for  ni  = 5d  and 

nH  ni  4f 
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5g_  14-16  self-explanatory  as  to  the  origin  of  the  values  we 

have  chosen.  With  these  and  the  <^r^  <Jr^  of  table  VI,  we  obtain  the 

final  normalized  factors  of  given  in  table  VIII.  These 

values  should  be  multiplied  by  the  corresponding  factor  in  table  VI, 
column  a,  for  each  N if  intensity  comparisons  between  ions  are  to  be 
made . 

TABLE  VII.  ENERGIES  IN  10^  cm’'  FOR  Af,  5d , AND  5g , DIFFERENCES 
l5d  = ^5d  - Ei,f  (AND  SIMILARLY  FOR  5g),  AND  RATIOS 
FOR  THE  TRIPLY-IONIZED  LANTHANIDES 


Ion 

C.  )'^ 

“5q^’5d 

Ce 

296.5 

246.8 

74.0 

49.7 

222.5 

4.48 

Pr 

31^. 

253.2 

76.0 

61  .2* 

238.4 

3.90 

Nd 

325.2 

254.8 

7b.  1. 

70.4 

248.8 

3.53 

Pm 

328.2 

256.6 

77.0 

71  .6 

25t  .2 

3.51 

Sm 

330.8 

258.3 

77.5 

72.5 

253.3 

3.49 

Eu 

3LI  .0 

260.0 

78.0 

81  .0 

263.0 

3.25 

Gd 

353.8 

261.5 

78.4 

92.3 

275.4 

2,98 

Tb 

318.6 

263.5 

79  0 

55.1 

239.6 

4.35 

Dy 

331  .9 

265.3 

79.6 

66.6 

252.3 

3.79 

Ho 

3LI  .6 

267.0 

80.1 

74.6 

261.5 

3.51 

Er 

392.6 

268.7 

80.6 

73.9 

262.0 

3.55 

Tm 

343.2 

270.5 

81  .2 

72.7 

262.0 

3.60 

Vb 

352.1 

272.2 

81.7 

79.9 

270.4 

3.38 

60  (1'>74), 


'■Sri 


^Measured  values.  See  footnote  a be2ow. 

L.  [Candor  Siuis  and  L.  J.  Nugent,  J.  Chem.  Phys 
1928,  table  II. 

^Calculated  by  subtracting  col.  d from  col.  a so  that  -£«; 

-^sf  - ^5d- 

^Calculated  by  - 0.3(~E^^) , which  gives  E^g(Nd)  • 

-76,400  cm~^--a  value  similar  to  that  obtained  from  the  Lu  III  spec- 
trum from  V.  Kaufman  and  J.  Sugar,  J.  Opt.  Soc.  Am.,  ^ (1971),  1693, 
as  follows:  From  * (^6s  ~ ^5g^  ~ =■  -169.049  and 

^6s  ~ ^5g  * -129.107  in  units  of  20^  cm~^,  then  nsg  - 4.9726.  Thus, 
16R/(nlg)^  . 71,000  cm-^'  -E^g 


cm 

77"^  ^ lanthanides . 

“K.  L.  V'anofer  Siuis  and  L.  J.  Nugent,  J.  Chem.  Phys.,  ^ (1974), 
1927,  table  I. 

^Calculated  by  subtracting  col.  c from  col.  a so  that  A5^  » ~^4f 


L.  Vander  Siuis  and  L.  J.  Nugent,  J.  Chem.  Phys.,  6£  (1  March 

1974),  1927;  J.  Opt.  Soc.  Am.,  64  (May  1974),  687. 

Kaufman  and  J.  Sugar,  J.  Opt.  Soc.  Am.,  6^  (1971),  1693. 

Sugar  and  J.  Reader,  J.  Chem.  Phys.,  59  (15  August  1973),  2083. 
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TABLE  VIM.  NORMALIZED  VALUES  FOR  TRI  PLY- I ON  I ZE  D LANTHANIDES 


N 

1 on 

a 

b 

c- 

d 

e 

1 

Ce 

2.01207 

9.50999 

0.0910822 

0.399768 

1 .975690 

2 

Pr 

1 .63399 

3.98929 

0. 0790669 

0.281959 

1 . 1 32910 

3 

Nd 

1 

2.89709 

0.0712982 

0,239639 

0.910776 

Pm 

1.39669 

2.72519 

0.0670561 

0.213609 

0.769900 

S 

Sm 

1.37931 

2.57925 

0.0636956 

0.192967 

0.661095 

6 

Eu 

1 . 23997 

2.21697 

0.0591970 

0.171118 

0.559329 

7 

Gd 

1 .08392 

1 .87091 

0.0599235 

0.151952 

0.975365 

8 

Tb 

1 .81988 

3.01990 

0.0616965 

0.163839 

0.991978 

9 

Dv 

1.50190 

2.91 173 

0.0575626 

0. 197127 

0.925229 

to 

Ho 

00 

0 

2.08196 

0.0598590 

0.1 35308 

0.377375 

) 1 

Er 

1 .35318 

2.03565 

0.0593990 

0.129668 

0.399897 

12 

Tm 

1.37552 

2.00798 

0.0591 681 

0.125325 

0. 327876 

13 

Yb 

1.25156 

1.77982 

0.0525539 

0.1181 70 

0.300956 

son: : 

Sornkil  1 z« 

■(/  by  using 

A5,-  (•  = d,  g) 

from  table  VII 

and  r ■■  r 

from 

CJblo  j 

•-Mii/ripJy  bii  corresfxjnill  Vciiues  in  t<ibia  H',  col.  d,  an.i  divide  by 
lOl'  to  obt<3in  actual  yaiut-s.  ; 

^^Normali  zcd  ' r ' r"*  j 

normal  izod  (r^‘ r j 

'^yor/vii  i zed  ■ r^' r ^ 

^noritial  i zed  r^'  r^?/l.i^q 
*^Hormalizivl  r>‘Cr'V.’-5^ 

4.  DISCUSSION  AND  CONCLUSIONS 

It  is  important  to  note  that  the  ostensibly  superior  method 
described  in  this  report  of  using  free-ion  wave  functions  to  determine 
phenomenological  for  Nd  in  CaW04  resulted  in  a fit 

(141  cm”^  compared  with  101  cm”l  rms)  between  phenomenological  and 
derived  inferior  to  that  obtained  previously^  (by  using 

Russell-Saunders  wave  functions  in  the  ESOH  method) . This  result  does 
not  indicate  that  the  discrepancy  in  the  fits  is  necessarily  due  to 

^R.  P.  Leavitt,  C.  A.  Morrison,  and  D.  E.  Wurtman,  Rare  Earth 
Ion-Host  Crystal  Interactions , 3.  Three- Parameter  Theory  of  Crystal 
Fields,  Harry  Diamond  Laboratories  TR-1673  (June  1975)  . 
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incorrect  phenomenological  B,  , but  that  further  theoretical  work 

;<m 

probably  is  required  to  determine  more  accurate  values  for  derived  B,  . 

km 

On  the  other  hand,  a comparison  of  the  Nd  results  with 

phenomenological  B,  determined  for  Er^"*"  in  CaWO., , also  using  free-ion 
km 

wave  functions,  resulted  in  an  approximate  constancy  with  respect  to  k 

and  m in  values  for  Jbj^^  (Er) /b^ (Nd)J  . Thus,  the  crucial  assumption 

in  the  theory  remains  so  far  intact — that  the  A,  are  impurity  ion 

km 

independent.  Without  this  assumption,  it  would  be  very  difficult  to 
continue  the  search  for  a theoretical  model  to  predict  impurity  ion 
spectra. 


There  are  several  areas  for  improvements  in  the  theory  that  readily 
come  to  mind.  First,  better  values  for  r and  may  be  determined  by 
including  electron  correlations^^  and  relativistic  corrections^®  in  the 
radial  wave  functions — effects  neglected  in  the  Hartree-Fock 
approximation.  Second,  one  may  determine  the  homogeneous  dielectric 
effect  of  placing  the  impurity  ion  in  a host  environment  and  also  the 
wave-function  overlap  effects  (which,  incidently,  might  in  part  be 
contained  implicitly  in  the  dielectric  calculation) . Third,  and 
probably  most  important,  the  dipolar  (and  possibly  higher  multipolar) 
contributions  to  the  A,  in  the  host  material  should  be  considered, 
since  their  contribution  has  already  been  shown  to  be  substantial. 


The  (table  II)  and  //S  values  (table  VIII)  , while  probably 

not  sufficient  for  predicting  accurate  spectra  and  electric  dipole 


C.  Morrison,  J.  Phys.  B:  Atom.  Molec.  Phys . , 6_  (November  1973), 

2205. 

B.  Lewis,  Relativistic  Calculations  of  and  Other 

^r  / Parameters  Encountered  in  Magnetic  Resonance  of  Rare-Earth  Ions  and 
Atoms,  Proc . XVIth  Congress  A ,M ,P .E .R .E . , Bucharest,  Rumania  (1  to  5 
September  1970) , (Publishing  House  of  Academy  of  Socialist  Republic  of 
Rumania,  1971)  . 

; A.  Morrison,  Solid  State  Commun,,  18  (1976),  153. 

t t — 


exclusively  in  the  following  10 
and  predict  spectra  for  the 


intensities,  nevertheless  will  be  used 
reports  in  this  series  to  interpret 
lanthanides  in  CaWOi,,  LiYF4,  YVOi,,  YPO4 , YASO4 , Y2SiBe2  07  , 
Y iAl 5O 1 2 (YAG) , Y3Ga50i2  (YGaG) , YAIO3 , and  LnP50j4,  where  Ln  is  any  of 
the  lanthanides.  These  calculations  will  enable  one  to  identify  the 
specific  areas  of  deficiency  in  the  present  theory  by  comparing  the 
predictions  with  the  mass  of  reported  data  in  the  literature. 
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APPENDIX  A. — HYDROGEN-LIKE  WAVE  FUNCTIONS  AND  EXPECTATION  VALUES 


From  equation  (16)  in  the  main  body  of  the  report  for  i j»  one 
obtains 


_ 2^1  /iL-V'5 


(A-1) 


5d 


W„ 


4f  3 \35j 

• Ir  (ff  “ r ♦ Ir  * “ S * 


5g  45 

From  the  same  equation  for  i / j , one  has 


(A- 3) 


C = 


y + X 


(A-4) 


1 - C 


l] 


(A-5) 


B = -C' 


I M)' 


(A-6) 


all  in  terms  of  the  exponential  factors  y and  x,  respectively,  for  the 
4f  and  both  the  5d  and  5g  wave  functions.  The  required  r^  expectation 
values  in  these  states  are  given  by 

. „2  (1-L^ 


(A-7) 


= W W 

\ ' 4f ,5d  4f  5d 


' "4f  (2y)9+P  ' 

r . . A !’■?  * B " P’- J . <A-8) 

L(y  ♦ (y  4 (y  4 
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4f,5g 


W W 
4f  5q 


(9  + p)  '■  (8  + p) 

, ,10+p  , ,9+p 

_(y  + x)  ^ (y  + x)  ^ 


(A-9) 


Thus,  all  of  the  normalizations  of  the  radial  wave  functions,  as 
well  as  the  coefficients  A,  B,  and  C defining  them  (according  to  eg 
(15a)  to  (15c)),  are  given  once  y and  x are  determined.  From  table  VI 
in  the  main  body  of  the  report,  one  has  for  Nd,  for  example. 


)^^(Nd)  ~ r‘*exp(-6r)  , 

(A-IO) 

^,(Nd)  ~ (r^  - 3r‘*  + 2r^)exp(-3r)  , 

bd 

(A-11) 

(Nd)  '•  (r^  - r‘*)exp(-3r)  . 

5g 

(A-12) 

In  retrospect,  it  might  have  been  better  to  choose  the  wave  function 
as  pure  hydrogenic  as  r^exp(-xr)  and  to  throw  the  deviation  from 


hydrogen- li)«e  behavior  in  the  wave  function,  because  the 
former — being  farther  removed  from  the  core — is  expected  to  be  more 
hydrogen-like . However,  since  the  wave  functions  are  used  to  calculate 
only  ratios  of  expectation  values  that  are  important  in  determining 
relative  electric  dipole  intensities,  the  ratios  are  not  expected  to  be 
sensitive  to  this  choice. 


A further  question  that  may  legitimately  be  asked  is  why  all  three 
radial  wave  functions  were  chosen  mutually  orthogonal.  It  is  true  that 
hydrogenic  radial  wave  functions  are  orthogonal  either  within  an  "n" 


shell^ij 

'nP'^n^'/ 

or  between  different  n 

shells 

if  S. 

and  i'  are 

equal. 

= 6 , . Thus , V is 

nn'  4f 

required 

to 

be 

orthogonal 

to 

and  not 

to  or  . However, 

bd  5g 

since 

the 

electron  orbital 

angular  momentum  K is  not  a good  quantum  number  so  close  to  the  core 
electrons,  an  orthogonality  between  the  approximating  radial  wave 
functions  is  probably  more  valid  than  if  such  a condition  were  not 
imposed . 
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